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Abstract 

In this paper we propose a class of exact solutions (stationary and non- 
stationary) of Einstein's field equations. We find that the space-time geometries 
of the solutions are non-vacuum and conformally fiat, whose energy-momentum 
tensors possess dark energy with negative pressure and the energy equation of 
state parameter w ~ —1/2. We also find that the time-like vector fields of the 
matter distributions of the solutions are expanding, shearing with acceleration 
and zero-twist. It is also found that, due to the negative pressure, the energy- 
momentum tensors violate the strong energy conditions leading to the repulsive 
gravitational fields of the space-time geometries. Energy-momentum tensors for 
the solutions also obey the energy conservation equations. From these physical 
properties of the matter distribution we may refer the space-times to as examples 
of exact solutions of the Einstein's field equations admitting dark energy with 
negative pressure. It is to note that the approximate sizes of the masses of the 
(stationary and non-stationary) solutions are less than (1/2) x 10~^° in Bousso's 
length scale r > 10^°. We also find that the surface gravities on the horizons are 
directly proportional to their respective masses. 

Keywords: Exact solutions; Einstein's equations; energy conditions; surface 
gravity; dark energy. 



1. Introduction 

Nowadays, it is found that the Wang-Wu mass function [1] plays a very important 
role in generating (embedded and non-embedded) exact solutions of Einstein's field 
equations. The mass function is expressed as a power series of the coordinate r as 

+ 00 

M{u,r)= ^ qn{u)r^, (1.1) 

n=—OD 

where qn{u) are arbitrary functions of retarded time coordinate u. The mass function is 
being utilized in generating non-rotating embedded Vaidya solution into other spaces 
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by choosing the function qn{u) corresponding to the number n [1]. Later on, the 
utihzation of these mass functions has been extended in rotating system and found 
the role of the number n in generating rotating embedded solutions of field equations 
[2]. The meaning of the power n in the expansion series (1.1) for the known spherically 
axisymmetric solutions are as follows [1, 2] 

(i) n = corresponds to the term containing mass of the vacuum Kerr family 
solutions such as Schwarzschild, Kerr; 

(ii) n = — 1 is equivalent to the charged term of Kerr family such as Reissner- 
Nordstrom, Kerr-Newman; 

(iii) n = 1 furnishes the term of the global monopole solution; 

(iv) n = 3 provides the de Sitter cosmological models, rotating and non-rotating. 

These values of n are conveniently used for stationary solutions (non-rotating and 
rotating). It is possible to obtain non-stationary Vaidya-Bonnor black holes, non- 
rotating [1] and rotating [2] when the values of n = and n = — 1 are summed up, 
since Vaidya-Bonnor metric describes a charged solution. When the above four values 
of n (= —1,0, 1,3) are together, one can find the charged Vaidya-de Sitter-monopole 
solution, non-rotating [1] and rotating [2]. Non-stationary rotating as well as non- 
rotating de Sitter cosmological models can also be obtained when n = 3 [3]. In fact 
the Wang-Wu power series expansion of the mass function turns out to be the most 
convenient method to generate new viable (embedded or non-embedded) solutions [2] 
of the field equations if one uses Newman- Penrose (NP) formalism [4]. From the above 
identification of power n we observe that the cases n = 2 and —2 in the mass function 
(1.1) have not been seen considered before so far in the scenario of exact solutions 
of Einstein's field equations. This is the main aim of the paper for generating viable 
solutions of physical interest and to analyze the nature of the matter distributions in 
the space-time geometry. Here we shall concentrate only the case n = 2. The other 
case n = — 2 may be discussed elsewhere. 

The important feature of the solutions with n = 2 proposed here is that the space- 
time metrics are non- asymptotically flat when r — >■ oo. However, they are flat at the 
origin r — > 0. Geometrically, it is acceptable that any curved space is locally flat. The 
idea of non-asymptotic flatness is in agreement with the Carter's suggestion [5,6] in 
generating exact solutions that "it is not necessary to assume asymptotic flatness nor 
make the assumption that there are no other Killing vectors than ^ = dt and rj = d^" 
[7]. Other examples of non-asymptotic space-times are (a) the de Sitter solution with 
cosmological constant and (b) rotating Kerr-NUT solution. 

The paper is organized as follows: Sections 2 and 3 deal with the derivation of 
a class of viable solutions (stationary and non-stationary) of Einstein's equations. 
We find that the masses of the solutions proposed here describe the gravitational 
fields of the space-time geometries and also determine the matter distributions with 
negative pressures, whose energy equations of state have the value —1/2. However, 
the energy-momentum tensors of the matter distributions with negative pressures do 
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not describe a perfect fluid, which can be seen in the next sections. These non- 
perfect fluid distributions are in agreement with the remark of Islam [7] - "it is not 
necessarily true that the field is that of a star made of perfect fluid" . We also find that 
each solution has a coordinate singularity with horizon. Consequently we discuss the 
areas, entropies and surface gravities at the horizons for the solutions. The existence 
of the horizons discussed here are also in accord with the cosmological horizon [8] of de 
Sitter space with constant A, which is usually considered to be a common candidate 
of dark energy with the equation of state parameter w = —1. The paper is concluded 
in Section 4 with reasonable remarks and evolution of the solutions with the physical 
interpretation. Thus we summarize the results of the paper in the following theorems: 

Theorem 1. An exact solution (stationary or non- stationary) admitting an energy- 
momentum tensor of a dark energy having negative pressure with equation of state 
parameter w = —1/2, is a non-vacuum, conformally fiat (NCF) space-time. 

Theorem 2. The energy-momentum tensor of the matter distribution in an NCF 
(stationary or non- stationary) space-time violates the strong energy condition leading 
to a repulsive gravitational force in the geometry. 

Theorem 3. The time-like vector fields of the matter distributions in the (stationary 
and non-stationary) NCF space-times are expanding, accelerating and shearing with 
zero-twist. 

Theorem 4. The surface gravity at the horizon of an NCF (stationary or non- 
stationary) space-time is directly proportional to the mass of the solution. 

Theorem 1 shows the physical interpretation of the solutions that all components 
of the Weyl tensors of the space-time metrics vanish indicating conformally flatness 
of the solutions. The energy-momentum tensors associated with the solutions admit 
dark energy having the negative pressures and the energy equation of state parameters 
w = —1/2. It is also found that the masses of the solutions not only describe the 
gravitational fields in the space-time geometries, but also measure the energy densities 
and the negative pressures in the energy-momentum tensors indicating the non-vacuum 
status of each solution. It is the assertion of General Relativity that "the space- 
time geometry is influenced by the matter distribution" [9]. Due to the negative 
pressure in the energy-momentum tensors, the violation of the strong energy condition 
is shown in Theorem 2 leading to a repulsive gravitational force in the space-time 
geometries. Theorem 3 shows the physical interpretation of time-like vector fields of 
the matter distributions. Theorem 4 indicates the existence of gravity on the horizons 
depending on their respective masses. It is to mention that Schwarzschild solution 
represents a vacuum, non-conformally flat space-time showing the difference from the 
non- vacuum, conformally flat (NCF) solutions discussed here. The presentation of the 
article is based on mathematical calculation for deriving exact solutions of Einstein's 
field equations. In this paper we utilize the differential form language developed by 
Mcintosh and Hickman [10] in Newman-Penrose (NP) spin coefficient formalism [4] in 
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(—2) signature as mathematical tool. 

2. An exact solution admitting non-perfect fluid 

We consider a line element of a general canonical metric in Eddington-Finkestein 
coordinate systems {u,r,9,(f)} 

ds^ = - '^M{u,r)^du^ + 2dudr - r'^{de'^ + si-n^e dcf)'^), (2.1) 

where M(m, r) is referred to as the mass function and related to the gravitational fields 
within a given range of radius r. Here u is the retarded time coordinate. 

By virtue of the Einstein's field equations Rab — {l/2)Rgab = —KTab associated 
with the above line element (2.1), we find an energy-momentum tensor (stress-energy 
tensor) describing the matter distribution in the gravitational field as 

Tab = f^Lh + '^p£{anb) + 2pm(am6), (2.2) 
where the quantities are found as 

H = -^M{u,r),u, p=^M{u,r),r, p = - ^M{u,r),rr (2.3) 
with the universal constant K = SttG/c^. Here ^a, Ua and nia are given as follows 

L = Sl na = ^{l-2rMiu,r)]5l + Sl 

ma = -^{5l+ismedt}, (2.4) 

where are real null vectors and rua is complex having its conjugate fha with the 

normalization conditions ia'n"' = 1 = —mafff' and other inner products are zero. From 
(2.3) we observe that there is no straightforward way for solving the non-linear Ein- 
stein's field equations with the mass function M{u, r) to generate a viable solution of 
physical interest. In order to have a meaningful physical interpretation of the solution 
(2.1) one has to consider some certain assumptions on the mass function M{u, r) as the 
line element having the energy-momentum tensor (2.2) with the quantities (2.3) has 
no reasonable interpretation of exact solutions. For instance, the Vaidya null radiating 
solution can be obtained when one assumes the mass function to be M(n, r) = M{u), 
leading to the condition p = p = in (2.3). Therefore, the mass function M(u,r) for 
obtaining viable solutions can, without loss of generality, be expressed in the powers 
of r as in (1.1) [1]. The above line element (2.1) with the mass function (1.1) includes 
most of the known solutions of Einstein's field equations, that can be seen with the 
identifications of power index n = —1,0,1,3, depending on the system (rotating or 
non-rotating) mentioned in the introduction above. For example, dc Sitter solution 
with cosmological constant A is obtained by setting qn{u) = A/6, when n = 3 and 
Qniu) = 0, when n 7^ 3 in (1.1), providing the mass function M{u,r) = (A/6)r^ and 
having energy density p* = A/K, and pressure p = —A/K in the non-rotating system 
[1]. Similarly, by assuming Qni^u) — when n — 0, and Qni^u^ — when n 0, 
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one can obtain the Schwarzschild solution with constant mass m in the non-rotating 
coordinate system. On the other hand, this choice of the power n = provides the 
Kerr vacuum solution in a rotating system. This shows the fact that, although there is 
no straightforward way of solving the field equations with the mass function M{u, r) 
in (2.1), the utilization of the Wang-Wu mass function (1.1) in the field equations, 
seems reasonable in generating viable solutions of physical interest. It is found that 
the case n = 2 in the power series expansion (1.1) has not been considered before in 
the scenario of exact solutions of Einstein's field equations. It is hoped that this case 
n = 2 may provide viable solutions of physical interest with reasonable interpretation 
of the matter distributions for both stationary as well as non-stationary space-times. 

Here we consider the case n = 2 in the mass function expansion (1.1) for generating 
a viable stationary solution of Einstein's field equations. Then we shall investigate 
the physical interpretation of the nature of the matter distribution in the space-time 
geometry. For this purpose we choose the Wang-Wu function qn{u) in (1.1) as 

, . ( m, when n = 2 , . 

= 1 0, when n ^ 2, ^^.5) 

such that the mass function takes the form 

+00 

M{u,r)= qn{u)r'' =mr'^, (2.6) 

n=— oo 

where m is constant and u = t — r is the retarded time coordinate. Using this mass 
function in (2.1) we find a stationary line element 

ds'^ = (1 - 2mr)du^ + 2dudr - r'^{de^ + sin^^^^^^ (2.7) 

where the constant m is regarded as the mass of a test particle present in the space- 
time and is non-zero for the existence of the matter distribution in the geometry. When 
u = constant, the surface is the future directed null cone. This line-element describes 

a stationary solution, and has a coordinate singularity at r = (2m)~^ describing a 
Lorentzian horizon. The line-clcmcnt (2.7) is certainly different from (a) Schwarzschild 
solution with g^u = 1 ~ 2M/r having singularity at r = 2M, and M being the 
Schwarzschild mass; and (b) de Sitter solution having a cosmological constant A with 
g^^ = 1 - (l/3)r^A singularities at r± = ±(3A)^/2. 

Now using the mass function (2.6) in (2.2) and (2.3) we find an energy-momentum 
tensor (stress-energy tensor) describing the matter distribution in the gravitational 
field of the space-time geometry (2.7) as 

Tab = 2p nt,) + 2p mf^Jh^) , (2.8) 



where the energy density p and the pressure p are found as 



The null vector Ha of (2.4) takes the form 

na = ^{l-2rm}Sl + Sl. (2.10) 

The equation (2.9) indicates that the contribution of the gravitational field to the 
energy-momentum tensor Tab having the negative pressure p is measured directly by 
the mass m of a test particle. Here from (2.7) and (2.9) we observe the key role of 
the mass m that it not only describes the curvature of the space-time in (2.7) but also 
distributes the matter field (2.8) present in the space-time geometry. Then we find 
the ratio of the pressure to the energy density as the equation of state for the solution 

This negative value of co is due to the negative pressure of the fluid. 

The energy-momentum tensor (2.8) obeys the energy conservation laws, given in 
the Appendix A in terms of NP spin coefficients: 

= 0, (2.12) 

which shows the fact that the metric of the line element is a solution of Einstein's field 
equations. The components of energy-momentum tensor may be written as: 

T: = T; = p, tI = tI = -p (2.13) 

for future use. We find the trace of the energy-momentum tensor T^b (2.8) as follows 

12 

T = 2{p-p) = —m. (2.14) 

Here it is found that p — p must always be greater than zero for the existence of the 
solution (2.7) with m 7^ 0, (if p = p implies that m will be vanished). It is to emphasize 
that the energy-momentum tensor (2.8) does not describe a perfect fluid, i.e. for a 
perfect fluid, one has T^^^^ = (p + p)uaUb — pQab with a unit time-like vector Ua and 
its trace T^p^^ = p — 3p, which is different from the one given in (2.14). 

Energy conditions: For the analysis of the energy conditions we shall introduce 
an orthonormal tetrad with a unit time-like vector ti" and three unit space-like vector 
fields v"-, w°', z°- using the null tetrad (2.4) with (2.10) such as 

Ua = -^{^a + no), Va = -^{ia-no), (2.15a) 

1 i 

Wa = -^{ma + ma), Za = - -j^{ma - nia) , (2.15b) 

with the normalization conditions Uau"" = 1, Vav"- = Waw"- = ZaZ°- = — 1 and other 
inner products being zero. Then the metric tensor can be written as 

9ab = UaUb - VaVb - WaWb - ZaZb- (2.16) 
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Now we consider a non-spacelike vector fields for an observer 

Ua = aUa + PVa + JWa + SZa, (2.17) 

where a, /3, 7 and S are arbitrary constants [11], subjected to the condition that 

f/"J7„ = - ^2 - 7^ - ^2 > 0. (2.18) 

Then the energy-momentum tensor (2.8) can be written in terms of the orthonormal 
tetrad vectors given in (2.15) as 

Tab = {P + P){UaUb - VaVb) - Vdab- (2.19) 

Now Tab W^U^ will represent the energy density as measured by the observer with the 
tangent vector U"" (2.17). Then we have the following energy conditions: 

(a) Weak energy condition: The energy-momentum tensor obeys the inequality 
TabU'^U^ > for any time-like vector [/" which implies that 

/9>0, p + p>0. (2.20) 



(b) Strong energy condition: The Ricci tensor for Tab satisfies the inequality 
RabVU^ > for all time-like vector C/«, i.e. TahU"-U^ > (l/2)r, which yields 

p>0, p + p>Q. (2.21) 

(c) Dominant energy condition: For all future directed vector V^^TabU^ should be 
a future directed non-space like vector. This condition is equivalent to 

p2 > 0, -p^> 0. (2.22) 



Here, we find that the value of p given in (2.9) does not satisfy the strong energy 
condition which has to be (2.21) as per the condition TabWU^ > (l/2)r with (2.19). 
This violation of the strong energy condition is due to the negative pressure, and may 
lead to a repulsive gravitational force of the matter field in the space-time. It is like 
the repulsive cosmological constant [8] to lead to the accelerated expansion of de Sitter 
space. This proves the stationary part of Theorem 2 above. 

We also find that the line element (2.7) of the space-time is conformally flat C"^^^ = 
0, i.e. all the tetrad components of Weyl tensor are vanished 

V'o = V'l = ^^2 = V's = ^^4 = 0. (2.23) 
The curvature invariant for the solution (2.7) is found as 



^abcd^ = — 




(2.24) 



which is regular on the 'singular' surface r = (2m) ^. This invariant diverges only at 
the origin. This indicates that the origin r = is a physical singularity. This shows 



7 



that the singularity of the solution (2.7) at r = {2m)^^ is caused due to the coordinate 
system, just like in Schwarzschild solution with the horizon r = 2m [12]. From (2.7), 
(2.9) and (2.23), we come to the conclusion of the proof of Theorem 1 stated in the 
case of stationary part above. 

Raychaudhuri equation: Here we shall analyze the nature of the time-like vector 
fields appeared in the energy- momentum tensor (2.19) for the stationary solution 
(2.7). The time-like vector is often considered to be the 4- velocity of a fluid. So 
the 4-velocity vector measures the kinematical properties of a fluid whether the fluid 
flow is expanding {Q = u":^ ^ 0), accelerating [iia = Ua-bu'' / 0), shearing aab 7^ or 
non-rotating {wab = 0) . We shall investigate the change of the volume expansion from 
the Raychaudhuri equation, such that we can understand how the negative pressure of 
the fluid affects the expansion of the solution. For this purpose we find the covariant 
derivative of the 4-velocity vector in terms of null tetrad vectors (2.4) with (2.10) 
as follows ^ ^ 

Ua;b = -^\m(ldb - nj,b) - -(1 + 2mr)m(„mj,)|, (2.25) 

where m is the mass of the solution. In deriving the above expression (2.25) we use 
the definition of u°' given in (2.15a). This expression of Ua-^h is convenient to calculate 
the (volume) expansion scalar 6 = u^.^ and acceleration vector iia = Ua-b'^^ as follows 

e = = (1 + 3r-m) (2.26a) 

Ua = ~m{ea-na) (2.26b) 
u% = -^(l-mr). (2.26c) 

We find that the vorticity tensor Wab = U[a;b] ~ ^[a^6] is vanished for the solution (2.7), 
and however, the shear tensor aab = '^{a;b) ~ '^{a'^^b) ~ {^/•^)&hab exists as 

(^ab = (44 + riaUb) - ^£{anb) + "^(a^b)} > (2-27) 

which is orthogonal to u"" {i.e., CFabU^ = 0). It is found that the mass m of the solution 
does not explicitly involve in the expression of aab but its involvement can be seen 
in the null vector Ua given in (2.4). However, the mass m directly determines the 
expansion as well as the acceleration iia as seen in (2.26). We find from (2.26) that 
the particle moving on the space-time geometry (2.7) follows the non-geodesic path 
of the time- like vector {ua-.bU^ 7^ 0). This establishes the key role of the expansion of 
the solution with acceleration. The vanishing of the vorticity tensor Wab = may be 
interpreted physically as saying that the matter field of the solution is twist-free (non- 
rotating) as mentioned earlier. This comes to the conclusion of the proof of Theorem 
3 for the stationary case. 

Now let us observe the consequence of the Raychaudhuri equation for the stationary 
solution (2.7). The Raychaudhuri equation is given by 

G = u\ + 2(«;2 - a^) - ^0^ + Rabu'^u'' (2.28) 
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where O = ©:au". The shear and vorticity magnitudes are 2(j^ = aah^^""^ and 2w'^ = 
Wabw'^^] Rab IS the Ricci tensor associated with the space-time metric gab- Then the 
Raychaudhuri equation for the twist-free time-Uke vector is found as follows 

e = -^(l + 2rm), (2.29) 

which takes the constant value = —2m? on the horizon r = {2m)^^, showing the 
constancy of the expansion along the time-like vector n". The impact of the negative 
pressure p given in (2.9) is taken care in Ricci tensor Rab in (2.28). The acceleration 
vector iia and its scalar ii":^ are negative, affecting the expansion rate in (2.28). It 
is to mention that the time-like vector is not a static Killing vector £ugab / 0, 
indicating the difference from any static time-like Killing vector {£^gab = 0) which 
has no expansion and shear. 



Surface Gravity: The line element of the stationary solution (2.7) has a horizon 
at r = (2m)^^. In this regard wc have to note the fact that the existence of the horizon 
is in accord with the cosmological horizon of de Sitter space with constant A [8,13], 
since it is regarded as a common example of dark energy with the equation of state 
parameter w = —1. Therefore, wc expect that it is highly important to observe the 
interpretation of the mass of the solution in connection with the area, entropy, surface 
gravity as well as the temperature for the horizon. So we find the area at the horizon 
r = (2m)~^ as 



7r /'2-7r 

-2 



^0 



A= / {gee g^)^ de d(l) =7rm-^ (2.30) 



and the entropy, from the entropy-area relation 5 = ^/4 [8], as 

S = ^Trm-^ (2.31) 

It indicates that the area and entropy will always exist for the solution with non-zero 
mass m. 

According to Carter [14] and York [15], the surface gravity k of a horizon is defined 
by the relation n^'V^n" = nn^, where the null vector n" in (2.10) above is parameterized 
by the coordinate u, such that d/du = n^Vb- Then the surface gravity is expressed in 
terms of NP spin-coefficient 7 as follows [16] 

K = nV6nX = -(7 + 7), (2.32) 

where 7 = — m/2. From this we find the surface gravity on the horizon r = {2m)~^ as 

K = m, (2.33) 

which shows that the surface gravity is directly measured by the mass, or in other 
words, it is directly proportional to the mass of the solution. This establishes the 
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proof of Theorem 4 in the case of stationary solution. Then the Bekeinstein-Hawking 
temperature for the model at the horizon is found as 



hK _ hm 
27rGkc 2^Gkc ^ ' ' 

where h is the reduced Planck constant, c the speed of light, k the Boltzmann constant, 
and G the gravitational constant. It indicates that the surface gravity and temperature 
of the horizon will never become zero for the existence of the stationary NCF solution 
(m 7^ 0). When the mass m = becomes zero, the hue element (2.7) will be a flat 
metric, the non-existence of the solution, and at this stage the surface gravity as well 
as the temperature will vanish. It is consistent with the property of flat space-time 
geometry, where there is no gravity, one cannot determine the surface gravity of the 
solution. It is noted that the surface gravity KSch of the Schwarzschild black hole is 
inversely proportional to the mass M as /csch = (4M)~^ on its horizon r = 2M [17]. 

Size of the mass: It is quite interesting to introduce a possible size of the mass of 
the solution discussed here. According to Bousso [13], stars are as distance as billions 
of light years, so r > 10^*^ and stars are as old as billions of years, t > 10^^. In 
this length scale, the size of the mass of the solution at the horizon r = (2m) ~^ may 
become 

m = ^r-^ < ^ X 10"^°, (2.35) 

which is slightly bigger than the size of the cosmological constant |A| < 3r^^ < 
3 X 10"-*^^° with the cosmological horizon ta = \/S/A [13]. The relation (2.35) may 
provide an example of a tiny test particle having a small size mass in the Universe. 

Kerr-Schild ansatz: The line element (2.7) can be expressed in Kerr-Schild 
ansatz on Minkowski flat background rjab as 

ds^ = dt^ - dr^ - r2(d6|2 + sm^9 d(f ) - 2mrdu^ (2.36) 

under the transformation t = n + r. This is the Kerr-Schild ansatz on the flat back- 
ground in spherical coordinate system 

gab = nab + 2Q£ah, (2-37) 

where Q = —mr, gab is the metric tensor of the solution and £a = is the null vector 
with respect to Qab and rjab given in (2.4). This confirms the fact that the stationary 
solution obtained here is an exact solution of Einstein's field equations. The line 
element (2.7) can also be written in the (t, r, 0, 0) coordinate system for future use as 
follows 

ds^ = (1 - 2mr)dt'^ - (1 - 2mr)-^dr'^ - r'^dO'^ - r^sin^^ (2.38) 



10 



under the transformation dt = du + {1/(1 — 2mr)}dr. This is a very famihar form of 
the hne element in General Relativity, and its determinant is \g\ = —r^ sin^ 9. In this 
coordinate system we can easily observe the singularity at the point r = {2m)~^. 

Physical interpretation of the matter distribution: In order to interpret 
the physical meaning of the negative pressure of the matter distribution in the NCF 
solution, we shall consider other space-times having negative pressure and the energy 
equation of state with minus sign. For instance, the cosmological constant A of the 
non-rotating de Sitter solution, whose energy-momentum tensor is T^^ = A.g'^f, is 
regarded as a common candidate of dark energy [13,19-27] having the negative pressure 
p = —A/K, and the energy density p = A/K with the equation of state parameter w = 
p/p = —1, where K = SttG/c^. T^^ violates the strong energy condition leading to 
the term - the repulsive (not attractive) cosmological constant A [8] . This equation of 
state w = — 1 is also satisfied for (a) the non-rotating non-stationary (time dependent) 
de Sitter with cosmological function A(u), whose Tab = —^'r'A(u)^uiah + A{u)gab with 
p = —A{u)/K, p = A{u)/K and u is the retarded time coordinate [3]; (b) the rotating 
de Sitter solution A{u) [3] having p = -r^A(it)(r^ -|- 2a? cos^ 9) / {K R^) and p = 
r'^A{u)/{KR^R^) at the poles 9 = 7r/2 or 37r/2, where R^ = r'^ + cos^ 9 with the 
non-zero rotational parameter a. It includes the case of the rotating stationary (time 
independent) de Sitter solution with constant A [2]. This indicates that in the study 
of dark energy problems one needs not concentrate only on the cosmological constant 
A, and that one can consider the non-constant cosmological function A{u) in the 
rotating as well as non-rotating de Sitter space-time geometries with the equation of 
state parameter w = —1 [2,3]. It is to mention the equation of state w for ordinary 
matter field distributions for better understanding the dark energy problem. That we 
have the equation of state for other ordinary matters having positive sign (i) w = 1 
for electromagnetic field having p = p = e^ / {K R!^ R?) of the Kerr-Newman black 
hole with the constant electrical charge e, and in Vaidya-Bonner radiating black hole 
p = p = e{uf' / {K R^ R?) with variable charge e{u) of retarded time coordinate u [2], 
(ii) w = 1/3 for radiation field with p—2p = [13,19-27]. This shows the fact that dark 
energy always has a minus sign in the value of the energy equation of state parameter; 
whereas the ordinary matter has a positive sign. The negative sign in the equation of 
state is an important property for any matter field distribution to be interpreted as a 
dark energy. 

The observations of luminosity-redshift relation for the type la supernovas [28-30] 
suggest that the missing energy should possess negative pressure p and the equa- 
tion of state w = p/p [31]. The negative pressure of the dark energy may be the 
cause of the acceleration of the present Universe. Although the dark energy has 
been sought in a wide range of physical phenomena depending on the value of the 
equation of state parameter w, (i) a quintessence field —1 < w < —1/3, (ii) the 
cosmological constant it; = — 1, (iii) a phantom field w < —1 [32,33], the nature of 
the dark energy still remains a complete mystery [33,34] without any proper space- 
time geometry, except the assumption of a line element of a perfect homogeneous and 
isotropic space-time having a compatible energy-momentum tensor of perfect fluid 
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T"^ = diag{p{t), —p{t), —p{t), —p{t)}, with p{t) and p{t) being the energy density and 
pressure of the matter distribution in the Friedmann-Robertson- Walker universe. 

Prom the above scenario of dark energy and with the findings here (i) the negative 
pressure (2.9), (ii) the energy equation of state with minus sign w = —1/2 (2.11), (iii) 
the violation of strong energy condition (2.21), and (iv) the accelerating expansion 
of the time-like vector (2.26) for the stationary NCF solution, we may regard the 
matter distribution (2.19) as an example of dark energy with negative pressure whose 
space-time geometry is the line element (2.7). 

3. Non-stationary solution with negative pressure 

In this section we shall develop a non-stationary (time dependent) version of the 
non-vacuum, comformally flat (NCF) space-time (2.7) discussed above. Since there 
is no straightforward way for solving the non-linear Einstein's field equations associ- 
ated with the mass function M{u,r) in (2.1) with (2.3) for a viable solution of non- 
stationary NCF space-time, we follow the Wang-Wu technique [1] as above. Therefore, 
we assume the Wang-Wu function having a variable mass m{u) as follows: 

+ 00 

M{u,r)= ^ qniu) r"" = m{u) r"^ (3.1) 

n=—OD 

when n = 2. Utilizing this mass function in the general canonical metric (2.1) we find 
a non-stationary line element as 

ds^ = {1 - 2rm{u)}du^ + 2dudr - r'^^dO'^ + sin^e^d^^). (3.2) 

where m{u) is considered to be a variable mass of a test particle in the non-stationary 
system. The above solution has a coordinate singularity at r = {2m{u)}'~^. Such 
a generation of a non-stationary solution from the stationary one (2.7) is acceptable 
in the framework of General Relativity that the non-stationary Vaidya null radiating 
solution with variable mass M{u) is a generalization of the stationary Schwarzschild 
vacuum solution with constant mass M. Similarly, the non-stationary de Sitter solu- 
tion [3] with a cosmological function A(u) having horizons at r± = ib{3A('u)}^/^ can 
be obtained from the stationary de Sitter model of constant A. 

Now using the mass function (3.1) in (2.3) we find the null density the energy 
density p and the pressure p as 

2 4 2 

associated with a energy-momentum tensor for a non-stationary matter distribution: 

Tab = I^Lh + {P+ P){UaUb - VaVb) - PQab- (3-4) 

Here la is the real null vector, Ua a unit time-like vector Uau"" = 1 and Va a unit 
space-like vector Vav"" = — 1 defined as in (2.15a). It is observed that the presence of 
the null density p in Tab is due to the non-constant mass m{u) in the field equations 
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showing the evolution of a non-stationary solution. From (3.3) it follows that the 
energy equation of state parameter w = p/ p has the value —1/2 with minus sign. If 
one sets the mass function m(u) to a constant m, the above energy-momentum tensor 
will become the one given in (2.19) for the stationary solution with the quantities p 
and p. The energy-momentum tensor obeys the energy conservation laws T'^^.f^ = 
as shown in Appendix below. It shows the fact that the non-stationary solution (3.2) 
is an exact solution of Einstein's field equations. It is also found that, due to the 
negative pressure, the Tab violates the strong-energy condition, 

^/^ + P>0, ^p + p + p>Q. (3.5) 

This violation indicates that the gravitational force of the non-stationary model is 
repulsive as in the case of stationary one (2.7) above which completes the proof of 
Theorem 2. However, Tab (3-4) with the negative pressure satisfies the week energy 
condition {i) p/2 + p > (ii) p/2 + p + p > and dominant energy condition (i) 
PP + p^ > 0, (ii) —pp + p^ — p^ > 0. Here we observe the difference between the strong 
energy conditions of the stationary (2.21) and that of the non-stationary (3.5) with 
the null density p. We also find that the non-stationary solution is conformally fiat 
Cl^ = 0, ^.e. 

■00 = V"! = '02 = '03 = ^"4 = 0. (3.6) 

This shows that the solution (3.2) with a mass function m{u) has the same char- 
acterization of conformally fiatness of stationary solution (2.21) with constant mass. 
Equations (3.3), (3.4) and (3.6) provide the proof of the non-stationary part of Theo- 
rem 1. It is also noted that the structure equation of the Riemann curvature invariant 
for variable mass (3.2) has a similar form of that of the solution (2.7) in (2.24) as 

160 

which is regular on the 'singular' surface r = {2m{u)}~^ . This invariant is divergent 
only at the origin r = 0, which is a physical singularity. 

We also have the expansion scalar and the acceleration vector iia with its scalar 
u"'.a for time-like vector Ua appeared in (3.4) 

@ = u% = -J^{l + Srm{u)} (3.8a) 
Ua = ~rn{u){ea - ria} (3.8b) 

u% = lrn{u),u - ^^^{1 - m{u)r}. (3.8c) 

However the shear tensor aab remains unchanged as in the stationary solution (2.24) 
and the rotation tensor Wab is vanished (zero-twist). This follows the proof of Theorem 
2 of the non-stationary case. Prom the Raychoudhuri equation, we have the rate of 
change of the expansion scalar as follows 

= lm{u),u - ^{1 + 2rm{u)}. (3.9) 



RabcdR'''''' = - — m\u), (3.7) 
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The surface gravity of the non-stationary solution at the horizon r = {2m{u)} ^ takes 
the form 

K = m{u) (3.10) 

which shows the proof of Theorem 4 of non-stationary part. We have also seen the 
evolution of the non-stationary solution with the mass function m{u) in (3.3), (3.8c) 
and (3.9). From (3.3), (3.5) and (3.8) we may regard the solution (3.2) as an example 
of a non-stationery space-time admitting an energy-momentum tensor of a dark energy 
with negative pressure having an equation of state parameter w = p/p = —1/2. 

4. Conclusion 

In this paper we develop a class of exact (stationary and non-stationary) solutions 
of Einstein's field equations describing non-vacuum and conformally flat space-times, 
whose energy-momentum tensors possess dark energy fluids with the negative pressure 
and the equation of state parameter w = —1/2. The most exotic property of these 
solutions is that the metrics describe both the background space-time structure and 
the dynamical aspects of the gravitational field in the form of the energy-momentum 
tensors. That is to say that the masses of the solutions play the role of both the 
curvature of the space-time (non-flat) as well as the source of the energy-momentum 
tensor with Tab 7^ (non-vacuum) measuring the energy density and the negative 
pressure. This indicates that when we set the masses of the solutions to be zero, the 
space-times will become the fiat Minkowski space with vacuum structure T^b = 0. 
In the case of Schwarzschild solution, the mass plays only the role of curvature of 
the space-time and cannot determine the energy-momentum tensor. That is why the 
Schwarzschild solution is a curved non-flat, vacuum space-time with T^b = 0. Here 
lies the advantage of the solutions (2.7) and (3.2) as non-flat and non-vacuum, space- 
times over the Schwarzschild. It is also to mention that the solutions discussed here 
provide examples of conformally flat space-times, while other examples of conformally 
flat solutions are the non-rotating de Sitter models with cosmological constant A [18], 
function A(u) [3] and the Robertson- Walker metric [12]. 

We find that the time-like vector of the source is expanding 7^ 0, accelerating 
iia (2.26b) as well as shearing aab / (2.27), but non-rotating Wab = 0. This means 
that the stationary observer of the solution does not follow the time- like geodesic path 
as Ua-bU^ / 0. Similarly, a non-stationary observer in (3.2) follows the non-geodesic 
path (3.8b). We also find that the energy-momentum tensors for the solutions violate 
the strong energy conditions. The violation of strong energy conditions is due to the 
negative pressure of the matter fields content in the space-time geometries, which can 
be seen in (2.21) and (3.5) above, and is not an assumption to obtain the solution (like 
other models mentioned in [33] ) . This violation indicates that the gravitational fields of 
the solutions are repulsive (as pointed out in [35]) leading to the accelerated expansions 
of the universe. The expansion of the space-time with acceleration is in agreement with 
the observational data [28-30] . It is also noted that the strong energy conditions for the 
stationary solution (2.21) associated with the stress-energy momentum tensor (2.19) 
and that of the non-stationary one (3.3) are different from that of the perfect fluid 
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(p > 0, p+3p > 0). This indicates that the strong energy condition is mainly depended 
upon the structure equation of a particular energy-momentum tensor. 

It is emphasized the fact that our approach in the development of the solutions 
here is necessarily based on the identification of the power n = 2 in the Wang-Wu mass 
function without any extra assumption. This identification of the power n = 2 in the 
mass function (1.1) has considered here for the first time, and not been seen discussed 
before in the scenario of exact solutions of Einstein's field equations. It is also noted 
that we do not consider the Friedman-Robertson- Walker metric, filled with perfect 
fluid, which is assumed to be the standard approach for the investigation of dark 
energy problem as mentioned earlier [13,19-26]. That is why the energy-momentum 
tensors associated with the solution (2.7) and (3.2) do not describe a perfect fluid. 
This fact can be observed from the trace T = 2{p — p) of the Tab for both stationary 
and non-stationary solutions. This non-perfect fluid distribution of the solutions is 
also in accord with Islam's suggestion that it is not necessarily true that stars are 
made of perfect fluid [7] . 

From the study of the above solutions, we find that the energy densities are only 
contributed from the masses of the matters. It is the fact that without the mass of 
the solutions, one cannot measure the energy density and the negative pressure of 
the energy-momentum tensors in order to obtain the energy equation of state w = 
— 1/2. This means that the negative pressures and the energy densities associated 
with the energy-momentum tensors (2.19) and (3.4) are measured by the masses that 
produce the gravitational field in the space-time geometries of the solutions. Hence, 
we may conclude that the stationary and non-stationary solutions may explain the 
essential part of Mach's principle - "The matter distribution influences the space-time 
geometry" [36]. It is emphasized that the equations of state parameters w = —1/2 
for the matter distributions (2.19) and (3.4) are belonged to the range — 1 < ■u; < 
focussed for the best fit with cosmological observations in [25] and references there in. 

The metrics appear singular when r = and r = (2m) ~^ for stationary and 
r = {2m{u)}~^ for non-stationary at a particular value of u. These values of r have 
special importance. At the origin r = 0, there is a physical singularity where the 
curvature invariants diverge as shown in (2.24) and (3.7); at r = {2m)~^ for stationary 
and r = {2m{u)}~^ for non-stationary, the invariants are well behaved and finite. 
Accordingly, we find areas, entropies as well as surface gravities at the horizons. It 
is found that the surface gravities given in (2.33) and (3.10) are directly proportional 
to their masses of the solutions. This indicates that the existence of the masses 
imply the existence of their surface gravities and the temperatures on the horizons. 
The existence of horizons is also in accord with the cosmological horizon of de Sitter 
space with constant A [8], which is considered to be a common candidate of dark 
energy with the parameter w = —1 [13, 19-26]. This parameter of equation of state 
is also true for both the cosmological constant A as well as the cosmological function 
A(u) of the rotating and non-rotating de Sitter solutions [3]. According to the length 
scale r > 10^*^ suggested by Bousso [13], we find the approximate sizes of the masses 
less than (1/2) x 10"^°, which are bigger than the size of the cosmological constant 
|A| < 3 X 10~^^° with the horizon va = v^3/A. It is noted that to the best of the 
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authors knowledge, the solutions are not been seen discussed before. We hope that the 
exact solutions (2.7) and (3.2) may provide examples of space-times admitting dark 
energy-momentum tensors (stationary and non-stationary) having negative pressures 
with the equation of state parameters w = —1/2 in the accelerated expanding space- 
time geometries. 



Acknowledgement 

The authors, Ibohal and Ishwarchandra acknowledge their appreciation for hospitality 

received from Inter-University Centre for Astronomy and Astrophysics (lUC AA) , Pune 
during their visit in preparing the paper. The work of Ibohal is supported by the 
University Grants Commission (UGC), New Delhi, File No. 31-87/2005 (SR). 



Appendix A: Energy conservation equations 

In this appendix we shall show the fact that the energy-momentum tensor (3.4) satisfies 
the energy conservation equations T°'''.jj = 0. These are four equations, which, using 
Newman-Penrose (NP) complex spin coefficients [4], can equivalently be expressed in 
three equations - two real and one complex. Hence, we find the following 

Dp = {p + p){p* + p*), (Al) 
Sp = pn* + {p + p){T-Tr), (A2) 
Dp + Vp = p{{p* + p*) - 2(e + e-)} ~ {p + p){p* + p*), (A3) 

where k* , p* , p* , r, vr, etc. are spin coefficients, and DjV and S arc the intrinsic 
derivative operators. These (A1-A3) are general equations for an energy-momentum 
tensor of the type (3.4). 

Now, in order to verify the conservation equations (A1-A3) for the components of 
T"^ with the quantities p, p, p given in (3.3) we present the NP spin coefficients for 
the non-stationary metric (3.2): 

K* = (7 = A = e = 7r = T = ^' = 0, 

P* = --, fi* = -^{l-2rm{u)}, (A4) 
r 2r 

^ = -a = cot 6*, 7 = -m(u). 

The intrinsic derivative operators are given as follows: 

D = r5„ = dr, 

V = n«a„, = du- ^{1 - 2rmiu)}dr, (A5) 



= m-da = ^\de + ^dA. 

■s/2r (. smd ^ ) 
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where ia, Ua and rua are the tetrad null vectors. The equations (Al) and (A2) are 
satisfied by using (A4) and (A5). By virtue of (3.3), (A4) and (A5), we find the left 
side of (A3) as 

4 2 

+ V/9 = -— m(n),u - -— r m{u){l - 2rm{u)}. (A6) 
Kr Kr'^ 

This can be shown equal to the right side of (A3) after using (3.3) and (A4). It leads 
to the conclusion of the verification that the energy-momentum tensor (3.4) satisfies 
the energy conservation equations T°'^.^ = 0. This indicates that the non-stationary 
line element (3.2) is an exact solution of Einstein's equations. It is also to mention 
that when m{u) sets to a constant m for the stationary solution (2.7), the energy- 
momentum tensor (2.8) satisfies the energy conservation equations (2.12). 
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